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Abstract. The accreting millisecond pulsar XTE J 1807-294 was observed as a Target of Opportunity (ToO) by XMM-Newton 
on March 22, 2003 after its discovery on February 21, 2003 by RXTE. The source was detected in its bright phase with an 
^ I observed average count rate of 33.3 cts s"' in the EPIC-pn camera in the 0.5-10 keV energy band (3.7mCrab). Using the 

earlier established best-fit orbital period of 40.0741±0.0005 minutes from the RXTE observations and considering a circular 
binary orbit as first approximation, we derived a value of 4.8+0.1 It-ms for the projected orbital radius of the binary system 
and an epoch of the orbital phase of MJD 52720.67415(16). The barycentric mean spin period of the pulsar was derived as 



o 

^ ■ 5.2459427±0.0000004 ms. The pulsar's spin-pulse profile showed a prominent (1.5 ms FWHM) pulse, with energy and orbital 

' phase dependence in the amplitude and shape. The measured pulsed fraction in four energy bands was found to be 3.1 ±0.2% 

(3 ! (0.5-3.0 keV), 5.4±0.4% (3.0-6.0 keV), 5.1+0.7 % (6.0-10.0 keV) and3.7±0.2% (0.5-10.0 keV), respectively. Studies of spin- 

profiles with orbital phase and energy showed significant increase in its pulsed fraction during the second observed orbit of the 
neutron star, gradually declining in the subsequent two orbits, which was associated with sudden but marginal increase in mass 
accretion. From our investigations of orbital parameters and estimation of other properties of this compact binary system, we 
Q ' conclude that XTE J 1807-294 is very likely a candidate for a millisecond radio pulsar 
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1 . Introduction transient, ultra compact binary systems having orbital periods 

U rp. .„. . , v^rc ncmoo/i a- < 80 minutes (Bildsten & Chakrabarty'MJT" Juett, Galloway 

The accreting mimsecond pulsar XTE J 1807-294 was discov- n i u . onm i a. . ^ onm, -ru u- 

tor & Chakrabarty 2003, Markwardt et al. 2002). The binary pa- 
- - ered as the fourth candidate of this class by 7^X77? (Markwardt, ^ . , rvrvirvrv^.* u-.j ca -.i- a 

o • ■ n o . ^^^o r, , ^, ^^^o rr,, . rameters impW siO.0 1 "0. 02 Mq whitc dwart donors With mod" 

Smith & Swank 20031 on February 21, 2003. The coher- » i . • u • i- • .u u- • n . i FwvtI 

■' erately high inclination in these binaries (Galloway et al. 120021 

ent pulsation of 5.245902ms was detected and subsequently n u . ^rvrv^ in i . i ^rvrv^ 

^ ^ Bildsten 2002; Markwardt et al. 2002 1. 
the orbital period of 40.0741 ±0.0005 minutes was established 

(Markward, Juda & Swank .2003) . confirming it to be the 

shortest orbital period out of five such pulsar systems dis 



Studies of the binary parameters of compact stars and their 
stellar companions provide a better understanding of the phys- 



, „ T,. ^1 ^ c -i^rrr- Tiorv-j ical conditions of this new class of objects, thek accretion 

covered so far. The spectral measurements of XTE J 1807- 

physics and possible evolutionary connections to the class of 



radio pulsars (Bhattacharya & van den Heuveri99r). With this 



294 showed a combination of absorbed black-body and a 

comptonization model with no absorption or emission lines . . , , , , ., , , ■ r , ^ , 

,„ ^ , r,,, , . ■„• J 1 in mind, we have done a detailed analysis of the ToO data of 

(Campana et al. 12003). The other accreting millisecond pul- ^„ , , . ^. , 

, f ^cAviiono /f a^<o /fo »r- a XTE J 1807-294 obtained by XMM-A^ewfon on March 22, 2003. 

sars known so far are SAXJ1808.4-3658 (2.49ms; Wijnands „. . . . . vt-t, ,1 cm oo^ u u a ■ 

o A r^i- rnsno vT-i- inc! one /t on iv ^ i The spin-period of XTE J 1 807-294 has been Confirmed, using 

& van der Klis 11998), XTE J1751-305 (2.30ms; Markwardt ,,.,.5,, , , , . , . . . ^, , 

& Swank^^, XTE J0929-314 (5.41 ms; RemiUard, Swank ^MM -Newton d^i^ and the projected semi-major axis of the bi- 

& Strohmayer2QQ2j, XTE J1807-294 (5.25ms; Markwardt, nary orbit derived usmg a technique appl^ 

„ „ „ , Rrron , ^, , . VT-T- Tioi,i TOO pulsars (Kirsch & Kendziorra 2003). We present new results on 

Smith & Swank I2003> and the most recent XTE J1814-338 i^,^,,^ . ,. , ^ ^ ,. , 

io 1- . 1 l AAA-M . T> . -1 J . J- £ . XTE J 1807-294 in this letter from our Studies of orbital param- 

(3.18ms; Sttohmayeret al. 2003V Detailed studies on the first ^ 

eters, orbital phase and energy dependent pulse profiles, which 



three candidates showed that these were short period. X-ray 



provide better insight into the physical conditions of this pulsar 
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2. XMM-Newton Observations 

XTE Jl 807-294 was observed by XMM-Newton on March 
22, 2003, starting at 2003-03-22 13:40:27 UT, with 
^obs-id 01579601 in revolution number 601. The EPIC-pn 
CCD camera was operated in Timing mode with the Thick 
Filter. For this analysis, data from a 9 293 s long exposure 
were used which covered almost four orbital periods of the 
source. 

The EPIC-pn camera provides high time resolution in 
Timing mode (0.03 ms), and Burst mode (7yus) with moder- 
ate energy resolution (E/dE = 10-50) in the 0. 1-15 keV energy 
band. This makes the pn-camera most suitable for simultaneous 
timing and spectral studies of millisecond pulsar sources. 

3. Data analysis 

3.1. Data reduction 

The data were processed with SAS 5.4.1. Event times were 
corrected to the solar system barycentre with the SAS tool 
barycen. In the Timing mode of the pn camera, a point source 
will be completely smeared out in the DETY direction. Hence, 
this mode provides only one-dimensional spatial information. 
Therefore, we used a nine column (37 arcsec) wide source 
extraction region centred on the source, consisting of CCD 
columns 33^1. CCD columns 03-1 1 were used for the back- 
ground extraction region. An average counting rate of 33.7 cts 
s"' was detected from the source in the 0.5-lOkeV band 
(3.7mCrab), while an average background of 0.4 cts s"' was 
detected. 

3.2. gfid search of the orbital parameters 

Using the best fit orbital period of 40.0741+0.0005 minutes 



(Markwardt, Juda & Swank I2003t . and assuming a circular or- 
bit, the relevant orbital parameters were determined through 
maximum epoch folding using 6 bin pulse profiles. To this 




Fig. 1. X grid search of the orbital parameters. See text for 
details. 

end, the photon arrival times with respect to the binary barycen- 



tre of all events were corrected for different values of projected 
orbital radius xq, and orbital phase to, where jcq was varied in 
seven steps of O.lOlight-ms and fo in 14 steps equivalent to 

1.0°, using fbinary con-ected = t - Xq sin(27r/Poi-b X (f - fo)). 

The initial parameter values were obtained by grouping 
events into 20 phase bins of the binary orbit and identifying 
two phases separated by 180° via ax^ epoch folding for each 
of the phase bins. 

The results of the two-dimensional parameter space search 
are shown in Figure The data were interpolated through 
Delaunay triangulation in order to better to visualise the;^'^ de- 
pendence on X() and fo. A maximum of 520 was found for xq 
of 4.8 +0.1 light-ms and fo of 52720.67415(16) (MJD). 

In order to test the assumption of a circular orbit, an analy- 
sis of possible orbital eccentricity of the neutron star orbit was 
carried out. However, no statistical significance for an eccentric 
orbit could be found and therefore the orbit is treated as circular 
in all further analysis. Table 1 summarizes relevant parameters 
of the accreting millisecond pulsar XTE J 1807-294 established 
so far 

Table 1. Measured parameters of XTE J 1807-294. 



Parameter 



Value 



Instrument 



RA 
DEC 

Orbital period 
Spin period 
Proj. orbital radius xq 
Orbital phase fo 



18:07:00.0 (12000) 
-29:24:30.0 (12000) 
40.0741(5) min 
5.2459427(4) ms 
4.8(1) light-ms 
52720.67415(16) (MID) 



Chandra' 
Chandra' 
RXTE' 
XMM-Newton^ 
XMM-Newton^ 
XMM-Newton^ 



'Markward, luda & Swank|2003 ^this work 



3.3. Pulse profile 

Using the best-fit orbital period of 40.0741+0.0005 min- 
utes (Markwardt, Juda & Swank 2003) and the orbital pa- 
rameters derived as described in Section 13.21 photon ar- 
rival times were corrected to the binary barycentre. We then 
derived the barycentric mean spin period of the pulsar of 
5.2459427+0.0000004 ms. By epoch folding the data with the 
derived period with respect to the epoch at MJD 52720.72457, 
we established spin pulse profiles in four different energy bands 
between 0.5-lOkeV as shown in Figure|2 These profiles show 
a very prominent pulse (1.5 ms FWHM). The shapes and rela- 
tive strengths of the profile varies with energy. We determined 
the pulsed fraction in these energy bands using the expression 



pulsed fraction 



|flux per bin - mean flux per bin| 
total flux 



The measured pulsed fraction was found to be 3.1+0.2% 
(0.5-3.0keV), 5.4+0.4% (3.0-6.0 keV), 5.1+0.7% (6.0- 
10.0 keV) and 3.7+0.2% (0.5-lO.OkeV), respectively, in the 
four different energy bands as in Figure|2l 
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Orbital Phase Bin 1 
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Orbital Phase Bin 3 
Pulsed Frac. = 0.05 1 ±0.004 



Orbital Phase Bin 4 
Pulsed Frac. = 0.035 ±0.005 



Orbital Phase Bin 5 
Pulsed Frac. = 0.034±0.004 



Spin Phase 



Fig. 2. Spin-pulse profiles of XTE J 1807-289 in different en- 
ergy bands. 

To study the variation of the spin-pulse profile with orbital 
phase, we grouped the data in six different phases covering the 
complete binary orbit, where phase starts at to - 44.08 s. The 
folded light curves coiTesponding to six orbital phases derived 
for the 0.5-lO.OkeV energy band are shown in Figure|3] The 
spin-pulse profile corresponding to orbital phase 3, shows a 
significantly higher pulsed fraction. We then performed a de- 
tailed analysis for each of the four observed orbits and detected 
a significantly higher pulsed fraction corresponding to orbital 
phase 3 of orbit 2, which gradually declines during orbits 3 & 
4. Figure l^gives the pulsed fraction calculated for six orbital 
phases for all the four binary orbits in the 0.5-3.0keV energy 
band. 

4. Discussion 

X-ray pulse profiles of XTE J 1807-294 showed both energy 
and orbital phase dependence during the XMM-Newton obser- 
vation. Pulse profiles of XTE J 1807-294 show orbital phase de- 
pendent changes and variation of pulsed fraction as quoted in 
Figure|3 These effects are also seen in other X-ray pulsars due 
to geometrical effects with respect to the line of sight and phys- 



Fig.3. Spin-pulse profiles of XTE J 1807-289 at different orbit 
phases in the energy range 0.5-lOkeV. In dotted red the or- 
bit averaged folded spin light curve from Figure|2]in the same 
energy range. 

ical conditions of the pulsar (Wang & Welter [11)811 1. The sud- 
den increase of the pulsed fraction value, as seen in Figure 01 
particularly during phase 3 in orbit 2, decreases gradually dur- 
ing subsequent orbits and is caused very likely by a sudden but 
marginal increase in mass accretion. This causes an increase in 
the pulsed fraction of the pulsar which then gradually decreases 
due to the reduction of accreting matter and hence the pulsed 
fraction when neutron star poles are visible at the line of sight 
corresponding to a particular orbit phase 3. Such changes in 
the pulse profiles and pulsed fraction due to increase in mass 
accretion rate during an outburst phase are commonly seen in 
X-ray binaries (Galloway et al. 120021 Mukerjee et al. [20001. 

Determination of binary orbital parameters enables us to 
estimate properties of the accreting binary system. The short 
orbital period of 40.0741 minutes of XTE J 1807-294 immedi- 
ately confirms that it is the most compact binary system out of 
the five accreting millisecond pulsars known so far. Estimations 
of the semi-major axis a sin / for this binary using our technique 
described above give 0.0048 It-s which enables us to estimate 
the mass function of this pulsar using the expression, 
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Fig. 4. Pulsed fraction values for 6 orbital phases indicated by 
different symbols for all the four orbits in the 0.5-3.0keV en- 
ergy band. 



/m = 1.074 X 10-^p2/^(fl sin;/ llt-s)^(Porb/ldr^Mo. 

The pulsar mass function /m - l-6x lO^^M© is the smallest 
presently known for any stellar binary. For a neutron star mass 
of 1 .4 Mq, the minimum companion mass would be 0.007 Mq 
for an inclination angle of / = 90°, suggesting that the com- 
panion of XTE J 1807-294 is in the regime of very low mass 
dwarf, of the order of 7 Jupiter masses. The nature of the com- 
panion in an ultracompact binary system depends on prior evo- 
lution scenario. Binary systems with Port < 80 minutes can 
be formed in at least two possible evolution as summarized by 
Deloye & Bildsten|2003|and the references therein. The first is, 
stable mass transfer on to a neutron star from an evolved main- 
sequence star or a He-burning star In this scenario, the main 
sequence star, when its core has nearly completed Hydrogen 
burning, is brought into Roche lobe contact due to loss of 
orbital angular momentum from magnetic braking which can 
evolve into a system with orbital periods comparable to those 
of accreting millisecond pulsars detected so far and can reach 
a minimum period of 10 minutes (Nelson, Rappaort & Joss 
UM^ Podsiadlowski et al lMHI Nelson & Rappaport [2003 1 
and when the system evolving through 40 minutes on the way 
out from the minimum period can have masses (0.01 Mq) more 
in line with the measurements (Nelson & RaDPaport l2003> . The 
second channel involves evolution through common-envelope 
phase during an unstable mass transfer from the white dwarf 
donor's progenitor star of either pure He or C/O and are most 
recently discussed in detail by Deloye & Bildsten (2003 1 . 

The composition of the donor in XTE J 1807-294 could be 
He/H, pure He or a C/O mixture depending on formation his- 
tory. Without further data, the composition of the donor can 
not be determined. However, one can examine the nature of the 
donor based on orbital information of the binary system by ap- 
plying models developed for ultra compact accreting millisec- 
ond pulsars by Deloye & Bildsten (2003 1. Their model con- 
siders low-mass < 0. 1 Mq White dwarf companion of arbitrary 
degeneracy of evolved He or C/O composition and central tem- 
perature of 10^-10^ as range of relevance for these objects. 



the corresponding central densities (10^ gcm^^ are such that 
Coulomb and thermal contributions to the equation of state pro- 
vide non-negligible corrections to degenerate electron pressure, 
affecting their M-R relations. The model thus based on equa- 
tion of state, adequately describe relevant physics and yields 
analytic description of the qualitative behavior of M-R rela- 
tions and how they are affected by Coulomb and thermal con- 
tributions to the equation of state. The model established a rela- 
tionship between the orbital inclination for ultra compact bina- 
ries and the donor's composition and central core temperature. 
The application of this model to three of the known accreting 
millisecond pulsar systems suggest that if the donors in all the 
three systems are He White dwarfs, the T=0 objects are allowed 
in XTE J0929-314 and XTE J1807-294 while XTE J1751-305 
requires a hot donor (Bildsten 2002 1. While C/O cold donor is 
possible for XTE J 1807-294, the other two systems required 
hot C/O donors. For XTE J1807-294 in particular, with mea- 
sured orbital parameters suggest that the companion for this 
ultra compact binary system has greater phase space available 
for a C/O donor than a He donor Because for XTE J 1807-294, 
the apriori probability greatly favors C/O donor in terms of al- 
lowed inclination as the probability that a donor is He is only 15 
% in this system (Deloye & Bildsten [2003> . Further evolution 
of this system will depend on its orbital inclination which in 
terms depends on donor mass, core temperature and its compo- 
sition. From known properties of this binary system, it is likely 
that at the end of its accretion phase, it could be a candidate for 
a millisecond radio pulsar 

Acknowledgements. The XMM-Newton project is an ESA Science 
Mission with instraments and contributions directly funded by ESA 
Member States and the USA (NASA). The German contribution of 
the XMM-Newton project is supported by the Bundesministerium fUr 
Bildung und Forschung/Deutsches Zentrum fiir Luft- und Raumfahrt. 

We thank an anonymous referee for his very constructive com- 
ments. 



References 

Bhattacharya, D. & van den Heuvel, E.P.J., 1991, Phys. Repts, 203,1 
Bildsten, L. & Chakrabarty, D., 2001, ApJ, 557, 292 
Bildsten, L., 2002, ApJ, 577, L27 

Campana, S., Ravasio, M., Israel, G. L., Mangano, V. & Belloni, T., 

2003, ApJ, 594, L39 
Deloye, C. I. & Bildsten, L., 2003, ApJ, 598, 1217 
Galloway, D. K., Chakrabarty, D., Morgan, E. H. & Remmilard R. A., 

2002, ApJ, 576, L137 
Juett, A. M., Galloway, D. K. & Chakrabarty, D., 2003, ApJ, 587, 754 
Kirsch, M.G.R & Kendziorra, E., 2003, ATEL 148 
Markwardt, C. B. & Swank, J. H., 2002, lAU Circ. 7867 
Markwardt, C. B., Swank, J. H., Strohmayer, T. E., In'T Zand, I. J. M. 

& Marshall, R E., 2002, ApJ, L21 
Markwardt, C. B., Smith, E. & Swank, J. H., 2003, lAU Circ. 8080 
Markwardt, C. B., Juda, M. & Swank, J. H., 2003, lAU Circ. 8095 
Mukerjee, K., Agrawal, P. C, Paul, B., Rao, A. R., Seetha, S., 

Kasturirangan, K. et al., 2000, A&A, 353, 239 
Nelson, L. A., Rappaport S. A. & Joss R C, 1986, ApJ 304, 231 
Nelson, L. A. & Rappaport, S., 2003, ApJ, 598, 431 
Podsiadlowski, Ph., Rappaport, S., & Pfahl E. D., 2002, Apl, 565,1 107 
Strohmayer,T. E., Markwardt C. B., Swank, J. H., et al., 2003, ApJ, 

596, L67 



M. Kirsch et al.: XTE J1807-294 



5 



Remillard, R.A., Swank, J. & Strohmayer, T., 2002, lAU Circ. 7893 
Wang, Y. M., Welter, G. L., 1981, A&A, 102, 97 
Wijnands, R. & van der Klis, M., 1998, Nature, 394, 344 



